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To all whom it may concern: 

Be it known Laurent Kieken, Matthew Neurock, Enrique Iglesia, and John Trenckle 
have invented certain new and useful improvements in 

METHOD AND SYSTEM FOR THE DEVELOPMENT OF MATERIALS 



m 



of which the following is a full, clear and exact description: 
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METHOD AND SYSTEM FOR THE DEVELOPMENT OF MATERIALS 

PRIORITY 



[0001] This application claims priority to, and is a continuation-in-part of, the U.S. patent 
application entitled A Knowledge-Based Process for the Development of Materials, filed 
August 6, 2001, having a serial number 09/921,929, the complete disclosure and drawings 
of which are hereby incorporated by reference. 

FIELD OF THE INVENTION 

[0002] The present invention relates generally to materials development methods, and, more 
particularly to a process for the rational development of materials used in chemical 
processes, including but not limited to heterogeneous catalysts. When applied to 
heterogeneous catalysts, this process may be referred to as a catalyst development engine 
(CDE). 



BACKGROUND OF THE INVENTION 

[0003] The still emerging, recent application of combinatorial chemistry to high-speed high 
throughput synthesis and screening of materials does not adequately address the 
commercial requirements that a catalyst must meet. Current combinatorial methods are 
based on random screening of large libraries of materials, prepared and evaluated under 
unrealistic conditions that are difficult to scale up. Thus, little useful knowledge is 
derived from such experiments to guide the selection of the next set of experiments or 
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materials and to scale up the material. A different approach for catalytic material 
discovery and development is needed in order to reduce the time to market which includes 
scalable high-throughput methods for catalyst synthesis and real-world conditions catalyst 
evaluation to accelerate generation of useful data coupled with a process that maximizes 
learning from these data and rapidly and efficiently identifies new material candidates. 
This knowledge driven process uses integrated scientific and empirical modeling tools to 
complement and mine experimental data in order to build predictive models that the 
scientist can use to guide material selection. This knowledge process for rational (as 
opposed to random) material discovery and development is the subject of this invention. 



use by the scientist to guide the selection of new materials and accelerate the rational 



development of materials. This system comprises a Knowledge Cycle (KC), a Testing 
Cycle™ (TC) and a knowledge management system. The KC comprises data-based and 
science-based modeling tools that are integrated in a knowledge management system, in 
order to maximize learning and enhance the scientist's decision-making capabilities for 
efficient experiment planning. 

[0005] Additional advantages of the invention will be set forth in part in the description 
which follows, and in part will be learned from the description, or may be learned by 




[0004] The present invention comprises a research process, preferably computer-assisted, for 



SUMMARY OF THE INVENTION 
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practice of the invention. The advantages of the invention will be realized and attained by 
means of the elements and combinations particularly pointed out in the appended claims. 

[0006] The above and other features and advantages are achieved through the use of a novel 
material development process as herein disclosed. In accordance with one embodiment of 
the present invention, a method of identifying chemical reaction mechanisms for a 
chemical process includes the step of specifying a reactant set that includes a plurality of 
chemical substances that may engage in a chemical reaction with one or more other 
substances in the reactant set. It also includes the step of specifying a plurality of possible 
q products that may result from the reaction of two or more of the substances included in the 

Q 

01 reactant set. It also includes the step of identifying a reaction mechanism set that includes 

m 

Hi a plurality of reaction mechanisms, wherein each reaction mechanism is a combination of 



4* two or more elementary steps representing the chemical process. It also includes the step 

Q of selecting a plurality of catalytic materials, where each catalytic material is associated 

rjj 

jjj with at least one of the reaction mechanisms in the reaction mechanism set, and each 

2 

ji j catalytic material is also associated with experimental data. It also includes the steps of 

associating a kinetic constant value with each elementary step of each reaction 
mechanism, as well as generating a kinetic model associated with each reaction 
mechanism and each catalytic material. Further, it includes the step of using a processing 
device to screen the reaction mechanism set by applying a goodness of fit test to the 
experimental data associated with each catalyst, eliminating the reaction mechanisms 
having a worst fit, and grouping the remaining reaction mechanisms associated with each 
catalytic material to provide a first reaction mechanism subset for each catalytic material. 
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[0007] Optionally, the method also includes the steps of selecting a performance variable and, 
for the reaction mechanisms contained in the first reaction mechanism subset, identifying 
one or more associated kinetic parameters to which the performance variable is most 
sensitive. 



[0008] The method of may also include the steps of using a processing device to calculate a 
modeled kinetic constant for a plurality of the elementary steps associated with a plurality 
of the reaction mechanisms using the processing device to screen the first reaction 
mechanism subset by eliminating the reaction mechanisms having associated kinetic 
constants that least closely relate to their corresponding modeled kinetic constants, and 
associating the remaining reaction mechanisms not eliminated in the second screening 
step with a second reaction mechanism subset. In this option, the calculating step may 
comprise using molecular modeling to calculate the modeled kinetic constant. Also, the 
option may include the additional steps of selecting a performance variable, and, for the 
reaction mechanisms contained in the second reaction mechanism subset, identifying one 
or more associated kinetic parameters to which the performance variable is most sensitive. 

[0009] In accordance with an alternate embodiment of the invention, a method of identifying 
materials for the performance of a chemical process includes the step of selecting a data 
set for a set of materials. The data set includes one or more dependent performance 
variables for a chemical process, as well as independent variables including, but not 
limited to, calculated or measured properties of the materials or preparation parameters 
relating to the materials. The method also includes the step of building a model that 
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correlates the dependent performance variables with one or more of the independent 
variables, as well as the step of identifying one or more of the independent variables 
having values that yield improved values of the dependent performance variables based on 
the results of the model built in the building step. Further, the method includes the step of 
identifying one or more new materials that are associated with the values of the one or 
more independent variables that yield improved values of the dependent variables. 



[0010] Optionally, in this embodiment, the step of building a model comprises the use of 
recursive partitioning. Also optionally, one or more dependent performance variables or 
one or more independent variables may comprise kinetic parameters that have been 
jjrj associated with reaction mechanisms in a reaction mechanism set. Also optionally, the 

m 

45 method may include the steps of applying a Monte Carlo kinetic simulation to calculate at 



a 



a 



least one modeled performance parameter for each material included in the material set, 
and selecting at least one materials class based on the results of the Monte Carlo 
simulation. 

[0011] Further, the method may include the steps of (i) selecting a selected reaction 
mechanism from a reaction mechanism set, wherein each reaction mechanism in the set 
comprises a combination of two or more elementary steps in a chemical process; (ii) 
applying a Monte Carlo kinetic simulation to calculate at least one modeled performance 
parameter for each material identified in the identifying step, wherein the simulation is 
associated with the selected reaction mechanism; and (iii) selecting at least one materials 
class based on the results of the Monte Carlo simulation. With this option, each reaction 
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mechanism in the reaction mechanism set may have been screened, using a goodness of fit 
test, to eliminate reaction mechanisms for which experimental data associated with 
reaction mechanism catalysts has been determined to have a poor fit. Each reaction 
mechanism in the reaction mechanism set may have been further screened to eliminate 
reaction mechanisms having associated kinetic catalysts that least closely relate to 
corresponding modeled kinetic constants. 

[0012] In accordance with an alternate embodiment, a process for the development of 
scalable, high-performance materials includes a computer-assisted knowledge cycle that 
uses at least one of (i) input from existing experimental data; (ii) correlations generated 
from at least one of experimental, theoretical, and/or modeling findings; and (iii) 
theoretical and modeling investigations to generate working hypotheses and suggested 
steps for at least one of experimental investigations and theoretical investigations to guide 
the search for better materials. 

[0013] Optionally, the knowledge cycle further may also include the use of kinetic modeling 
to guide catalyst development. The knowledge cycle may also include the use of machine 
learning methods to guide catalyst development, as well as using kinetic Monte-Carlo 
simulation to screen catalytic surfaces for catalytic performance. 

[0014] There have thus been outlined the more important features of the invention in order 
that the detailed description that follows may be better understood, and in order that the 
present contribution to the art may be better appreciated. There are, of course, additional 
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features of the invention that will be described below and which will form the subject 
matter of the claims appended hereto. 

[0015] In this respect, before explaining at least one embodiment of the invention in detail, it 
is to be understood that the invention is not limited in its application to the details of 
construction and to the arrangements of the components set forth in the description or 
illustrated in the drawings. The invention is capable of other embodiments and of being 
practiced and carried out in various ways. Also, it is to be understood that the phraseology 

p and terminology employed herein, as well as the abstract, are for the purpose of 

□ 

rj description and should not be regarded as limiting. 

81 

DO 
FU 

%[ [0016] As such, those skilled in the art will appreciate that the conception upon which this 

disclosure is based may readily be used as a basis for designing other structures, methods, 
and systems for carrying out the several purposes of the present invention. It is important, 
therefore, that the claims be regarded as including such equivalent constructions insofar as 
they do not depart from the spirit and scope of the present invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



[0017] The accompanying drawings, which are incorporated in and constitute a part of this 
specification, illustrate one embodiment of the invention and together with the description, 
serve to explain the principles of the invention. In the drawings: 

[0018] FIG. 1 is a schematic diagram showing a preferred embodiment of the catalyst 
development engine. 

[0019] FIG. 2 is a schematic diagram showing the underlying architecture of the catalyst 
development engine. 



[0020] FIGs. 3a and 3b illustrate exemplary and optional steps included in preferred 



[0021] FIGs. 4a through 4c illustrate exemplary results of microkinetic analysis for NO 
decomposition on Pt/Al 2 0 3 catalysts using the present invention. 

[0022] FIG. 5 is a plot of the apparent rate constant for NO decomposition versus oxygen 
affinity for various catalysts at two temperatures, in accordance with the example started 
in FIGs. 4a through 4c. 

[0023] FIG. 6 illustrates exemplary data for the identification of key catalyst properties. 
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embodiments of a knowledge cycle. 
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[0024] FIG. 7 is a plot of the apparent rate constant for NO decomposition versus the standard 
enthalpy of formation of oxides, in accordance with the example started in FIGs. 4a 
through 4c. 

[0025] FIG. 8 is a list of elementary steps for NO decomposition and oxidation, in accordance 
with the example started in FIGs. 4a through 4c. 

[0026] FIG. 9 is an energy diagram for the elementary steps listed in FIG. 8. 

M 

Q 
D 

^ [0027] FIG. 10 is a plot of the calculated binding energy of adsorbed oxygen on Pt(100) as a 

00 

fjj function of coverage, in accordance with the example started in FIGs. 4a through 4c. 
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[0028] FIG. 1 1 is a plot of the calculated turnover frequencies of N 2 and N0 2 formation on 
Pt(100), in accordance with the example started in FIGs. 4a through 4c. 

[0029] FIG. 12 is a plot of the calculated turnover frequencies of N 2 and NO, formation on 
Rh(100), in accordance with the example started in FIGs. 4a through 4c. 

[0030] FIG. 13 shows DFT calculated hgfcts of adsorption for O*, N*, and NO* and the NO 
activation barrier over several different transition metal surfaces, in accordance with the 
example started in FIGs. 4a/mrough 4c. 
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[0031] FIG. 14 represents the variations of the calculated NO activation energy as a function 
of the structure of ideal Pt single crystal surfaces, in accordance with the example started 
in FIGs. 4a through 4c. 

[0032] FIG. 15 shows the DFT calculated transition state for the activation of NO over 
Pt(100), in accordance with the example started in FIGs. 4a through 4c. 

[0033] FIG. 16 illustrates an exemplary computer suitable for carrying out the automated 
functions of the present invention. 

[0034] FIG. 17 illustrates the internal hardware of the exemplary computer of FIG. 16. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0035] The present invention is a method of discovering, developing, and optimizing 
materials. The method may be applied to the development of new heterogeneous catalytic 
systems. Reference will now be made in detail to the present preferred embodiments of 
the invention, examples of which are illustrated in the accompanying drawings. 

[0036] This section provides an overview of/me catalyst development engine (CDE) of the 
present invention and discusses the tools and methods required for data analysis and 
mining, predictive model development and catalytic materials selection. The topics 
discussed under this headjrfg are the key components of the novel CDE depicted in FIG. 1. 
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The engine provides a rapid approach to the rational development of scalable 
heterogeneous catalysts and of lygn-performance solid materials. The CDE includes three 
main components that are^mtimately connected in an integration framework: (1) a testing 
cycle (TC) 20 for experimental data generation; (2) a knowledge cycle (KC) 22 for data 
analysis and^pubsequent experiment planning; and (3) a knowledge repository 24 (e.g., a 
database^ 



[0037] The primary function of the TC 20 is to rapidly generate^informative data for the 
experimental evaluation of new catalytic materials and/also generate data on well- 
characterized systems and for probe reactions of materials properties such as, but not 
restricted to, acidity, basicity, reducibility, etc. wr use in the KC 22. The primary 
functions of the KC 22 are (1) to generate working hypotheses relating performance to 
key catalyst properties that will guide the/earch for even better catalysts and (2) to 
generate fundamental structure-property ^relationships that will be instrumental in the 
selection of catalysts for further experimental evaluation in the TC 20. Both cycles run 
concurrently and feature theoretical And experimental activities that are highly integrated. 
The knowledge repository expands after each iteration of the CDE. The knowledge 
repository inventories materials properties and catalyst performance data for a 
comprehensive database of ^materials and for characteristic probe reactions of materials 
properties. It contains yooth experimental and theoretical data. Such an approach 
accelerates developme/t and scale-up of new materials without the impediments 
introduced by conventional combinatorial approaches based on randomly selected 
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materials. Preferred components apdactivities of the CDE are described in more detail in 
the following table (Table l),^nd illustrated in FIG. 2. 



Table 1. Catalyst Development Engine - Research Process Description 


Statement of ProbU 
reasonably and adec 


?w: Entails a concise statement clearly defining project goal(s) in terms that can be 
uately addressed with the tools available. 


Prior Data 


All pertinent available data, independent of source, are collected, evaluated and stored 
in the knowledge repository in a format that is usable in the CDE. 


Testing Cycle: In this cycle, the emphasis is on a rapid selection, preparation, and evaluation of potential 
catalysts. 


Catalyst & Test 
Selection 


The catalyst target is selected for evaluation, which includes determining what 
performance tests are meaningful, e.g. rate, selectivities, conversions, heats of 
adsorption, etc. 


P mp ri mpnfnl 

Design 


Thp c« a lf»r , t* i H pfltnlvctc nr*» nrpnnrpH and fvnlimtpH arrnrHino tr\ c\r\ pvnpri mpntal H*»cion 
l nc aeieeieu uuiuij'iio die jjiejjuieu aim cvaiuai&u avuvjiuing tu an CApci iniciiiai ucMgn 

matrix. This ensures that the significant variables are effectively examined with the 
minimum number of experiments, while maximizing the information gathered. 


Synthesis, and 
Evaluation 


The catalysts are prepared using validated methodologies and HTP tools. They are 
synthesized and evaluated for key properties and performance as prescribed by the 
experimental design. 


Knowledge Cycle: 1 
(1) Discriminate bel 
or empirical modeli 


"his cycle uses literature, Testing Cycle data, theoretical chemistry and simulation data to 
ween reaction mechanisms; (2) Identify key measurable catalyst properties via molecular 
ng; and (3) Derive structure-property models to guide catalyst search. 


Screen Reaction 
Mechanisms 


• Microkinetic modeling is applied to identify critical mechanistic steps that limit 

nprfnrmanfp nnri nccrw^iatpH Vinptif* naramptprc Ttc fiinptionc inplnHf** 
pci iKji mauce aiiu assu^iaitu A.uieiie paiaiueiei o. lis luiieuvjua nieiuue. 

Generate mechanisms based on reactant, product, and possible intermediates 
Discriminate between reaction mechanisms based on goodness of fit with 
experimental data and with realistic constraints on kinetic parameters based on 
measurements or molecular modeling calculations. 

Perform sensitivity analysis on kinetic parameters to identify important ones. 

• Mechanisms that can reproduce general trends in existing data are studied in more 
detail with molecular modeling to further evaluate their validity. 


Identify Key 
Catalyst 
Properties 


• Experimentaltheoretical and simulation data are modeled with empirical machine 
learning algorithms to generate hypotheses relating catalyst descriptors including 
properties (e.g. binding energies, activation energies, other surface properties) and 
synthetic parameters to key kinetic parameters or to performance 


Guide Selection 
of Material 

Candidates for 
Subsequent 
Experiments 


• Electronic structure, composition/structure, particle size, support, and promoters 
can affect the surface properties of catalysts. Developing trends that relate catalyst 
descriptors to these properties will assist the scientist in deciding what "knob" to 
turn in order to tune the catalyst properties. 

• Screen virtual surfaces for performance using Monte-Carlo simulation using semi- 
empirical and ab initio computational methods for kinetic parameter estimation. 
Develop guiding trends across materials classes. 

• Develop optimal lead catalytic systems using predictive machine learning methods 
such as recursive partitioning, neural networks or genetic algorithm to estimate 
performance of new materials from catalyst descriptors 


Iterations: Using the knowledge gained from the current cycle a determination is made on how to proceed; 
either by preparing another round of targets or by terminating this phase of the project because the catalyst 
meets the project goals or it is a dead end and the project continuation is evaluated. 


Solution to Problem: When a new catalyst that meets the performance criteria is identified, it is considered 
for scale-up and commercial development 
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[0038] The process outlined in Table 1 is a nover decision-support tool for the researcher to 
use interactively in the development of new catalytic materials. The CDE integrates 
various modeling tools - computational chemistry tools, kinetic modeling tools, and 
machine learning tools - wittr a knowledge management system to house the significant 
experimental and theoretical knowledge for its efficient future use in the guided 
development of new^atalytic materials. 

[0039] FIGs. 3a and 3b provide additional detail about the knowledge cycle element of the 
present invention. Referring to FIG. 3a, the performance target or targets and desired 
reaction(s) are identified (step 70). Theoretical and experimental data are also gathered 
(step 72), from existing sources, from new research, and/or from the knowledge cycle 
process itself. The knowledge cycle can use either a microkinetic modeling approach 
(step 74) or a machine learning approach (step 80). Optionally, a combination of the two 
approaches may be used. 

[0040] With the microkinetic modeling option (step 74)/fne cycle includes the specification 
of reactants, products, and intermediates (step 7£F), the generation of reaction mechanisms 
(step 76), screening of mechanisms (step/77), and the identification of critical kinetic 
parameters (step 78). The screening step (step 77) is further illustrated in FIG. 3b. 
Referring to FIG. 3b, based onya microkinetic model 90 and experimental data 92, a 
goodness of fit test is perfonned (step 92) between the experimental data and the data 
generated by the model to Eliminate the reaction mechanisms for which the fit is not good. 

13 
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A second cut is taken by using kinetic paramete^calculated with molecular modeling 
(step 97), and assessing the reasonableness/of kinetic parameters associated with the 
mechanisms (step 96). The relative sensitivity of the performance of the catalyst to the 
kinetic parameters is also assessed (step 98), with a preference expressed for the 
parameters having the highest levels of sensitivity. 

[0041] Returning to FIG. 3a, in addition to the microkirf£tic modeling, databases of 
information gathered can be mined using one or more^of many types of machine learning 
and pattern recognition techniques to identify correlations between performance, 
properties, and preparation variables (step 8(f). These correlations become the basis for 
the selection of materials. These potential materials can be screened further using Monte 
Carlo kinetic simulation (step 88). In this step, catalytic performance of a virtual catalyst 
candidate can be calculated. Usmg one or a combination of these options, materials are 
selected for synthesis and testing (step 84), experiments are designed (step 84), and 
performance and properties are measured in a testing cycle (step 86). 

[0042] The following sections describe, in additional detail, the modeling tools that are 
suitable for treating heterogeneous catalyst problems: the computational chemistry 
methods, the microkinetic and Monte-Carlo modeling techniques for performance 
prediction, and the machine learning algorithms that together make the CDE™ approach 
novel and viable. 

Computational (Modeling) Methods 
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[0043] Current state-of-the-art molecular modeling of catalysts/involves either ab initio 
quantum chemical methods or semi -empirical method?: Ab initio quantum chemical 
methods are first-principle solutions to the Schrodkfger equation, which make only a few 
very basic approximations. Because of the/mgh level of detail and accuracy, they are 
applicable to essentially any materia^ system but they require high computational 
intensity. They provide the best accuracy for a wide range of chemical and catalytic 
properties. These methods are' used to examine the chemisorption and reactivity of all 
reactants, intermediates, ana products at different surface coverages in order confirm the 
reaction mechanism and provide the parameters necessary for detailed microkinetic 
models (adsorptior/energies, activation energies, rate constants). In addition, a large set 
of chemical anja catalytic properties for different catalytic systems are evaluated in order 
to develop/a set of property libraries that are used to establish correlations (metal- 
adsorbatp bond strengths, site acidity, site basicity, oxidation potential, etc.) and 
hypotheses tests. 



[0044] Semi-empirical methods include both empirical theoretical models, such as Bond 
Order Conservation, and semi -empirical quantum chemical methods. Semi -empirical 
quantum chemical methods attempt to minimize the computational efforts of solving the 
Schrodinger equation by approximating computationally limiting multi-center integrals. 
These approaches range from.very approximate calculations such as extended Huckel 
calculations, which simply ignore the contributions from multi-center integrals, to very 
approximate methods using low computational times on the order of a few minutes. More 
sophisticated methods treat the multi-center interaction integrals using empirical 
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information. They typipally use experimental data or ab initio results in order to 
parameterize the inj^grals. 

[0045] In terms of ab initio methods, Density Functional Theory (QPT) calculations provide 
the state-of-the-art for modeling catalytic systems and J#ey can be used to predict 
chemisorption energies, overall reaction energies, attivation barriers, and chemical 
properties descriptors which target key catalytic properties, such as acidity, basicity, and 
reducibility. Gradient-corrected DFT method^ have proven to be the most robust and 
accurate methods for the prediction of large/transition -metal systems due to their ability to 
explicitly treat electron correlation. Th/active surface sites can be approximated by using 
either a metal cluster model to represent the coordination and bonding of the active metal 
particles, or a periodic slab mod£l that represents larger exposed surfaces. The accuracy 
of DFT methods in predicty^g energetics for transition metal systems is +/- 5 kcal/mol. 
Although quantitative accuracy is important, the ability to predict the relative trends 
across the periodic ta/le is more critical for the identification of useful catalyst targets. 
This enables scannAg of a range of unknown systems for optimal metal-adsorbate bond 
strengths and surface reactivity. 



[0046] These detailed calculations can be used tgf understand and predict the trends in both 
chemisorption and surface reactivity apfoss the periodic table. Non-local cluster and 
periodic slab DFT calculations an? used to compute the binding energies of maleic 
anhydride on the 111 surface/of Pd, Re, Au, Pd/Re, Pt, Pd/Mo, and Pd/Au. Basic 
^ concepts from frontier mdecular orbital theory are subsequently used to construct 
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general model that is able to predict a priori the outcome^? the more detailed ab initio 
calculations. These models provide an understanding ofrthe controlling factors that govern 
chemisorption as well as a knowledge-driven approach to the screening of multi-metallic 
and metal oxide systems by simply computing particular chemical descriptor (in this 
case the center of the d-band at the surface/layer of the metal). This model is further 
extended to predict activation barriers yfor different elementary steps in the overall 
catalytic path. As such, studies show/that the reactivity of the surface scaled with the 
relative location of the d-band center of the metal with respect to the Fermi level. Also, 
the reaction kinetics over different (111) surfaces is found to correlate quite well with 
changes in the center of th/ d-band. There is an optimal trade-off between different 
reactions. As the metal-^adsorbate bond becomes too weak it becomes difficult to 
dissociate hydrogen on jne surface. This shuts down the overall reaction. 




[0047] While DFT quantum mechanical method arp'appropriate for establishing the reaction 
mechanism, and for accurately calculating^roperties for a large number of systems, the 
calculations can still require 1-2 days per simulation. In order to offer daily guidance and 
be able to lead high-throughput screening efforts, a much faster method needs to be 
developed. The present invention uses a variety of semi -empirical methods that run for a 
few minutes to determine specific catalytic properties within reasonable accuracy. These 
methods, driven by the scientist, can be coupled with advanced machine learning methods, 
that are described later, in order to develop a potential set of lead catalytic systems which 
meet the target properties. 
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[0048] Another semi-empirical method is the Atomic Superposition Electron Derealization 
Molecular Orbital Theory (ASED -MO): ASED may be used to describe many qualitative 
features and, in some cases, evepfquantitative assessments in catalysis. The applicability 
and accuracy of the approapn depends on developing appropriate parameters. 



[0049] The present invention generates very detailed m initio libraries as well as empirical 
databanks, which can be used to parameterize/mese systems more accurately. By way of 
simple comparison, we examined how wal we could do in terms of producing trends by 
comparing simple ASED calculatioj/which were optimized to detailed DFT calculations 
and found that the results showlne correct trends to within 30% for NO, O, and N binding 
energies over a set of different metals. The application of structural optimization as well 
as refined regression ofrme parameters should improve the accuracy to 10-20%. 



[0050] In addition to refined ASED models, we/nave developed a more approximate Bond- 
Order Conservation (BOC) method. Whnc BOC has some ties back to quantum chemical 
descriptions of bonding, it is a very simple approach that does not attempt to solve the 
Schrodinger equation. The approach offers some predictions to reactivity. The present 
inventors have found that/some of the known shortcomings of this method can be 
overcome by determining interaction parameters from first-principle DFT methods. Also, 
tight binding or singl/sCF DFT methods are used and evaluated for the most appropriate 
compromise between accuracy and speed. 



Microkinetic Analysis 
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[0051] Microkinetic modeling incorporates the basic surface chemistry, i.e., elementary steps, 
in the kinetic description of a catalytic reaction. Such a/kinetic model is a very useful tool 
to compare and extrapolate the performance of different catalytic materials at various 
process conditions. It is used to eliminate postulated reaction mechanisms that are not 
consistent with experimental data. Giverf a plausible mechanism, microkinetic analysis 
can also point to the often few critical/kinetic parameters that are sufficient to describe the 
overall rate of the catalytic process. This is accomplished by performing a sensitivity 
analysis of catalytic performance on all kinetic parameters. Knowledge of the critical 
kinetic parameters and associated steps can provide the scientist with insights on how to 
modify materials in ovd/fr to increase the overall performance. 



[0052] In the Catalyst Development Cycle, the microkinetic analysis may^nclude one or more 
of the following steps: 



(1) Mechanism enumeration: Given a catalytic reaction, various possible 

mechanisms are enumerated based on the list of observed products and 
possible intermediates, list of pfausible elementary steps, experimental 
kinetic data on a training se/of materials, literature data, and chemistry 
rules. For example, steps/that require the simultaneous reaction of more 
than three reacting sn6cies or the breaking and forming of too many 
chemical bonds are /sually not elementary and can be excluded. Reaction 
pathways can al^o be computer-generated by assuming again that an 
elementary ste/involves only a small number of changes in the bonding of 
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the reactants, i.e., three or four changes in the connectivity of the re^ctants 
by bond cleavage and formation. 
(2) Mechanism discrimination: There is usually a deficit of information for 

the values of kinetic rate constants of elementary steps, /specially for the 
activation barriers. Good estimates of pre-exponential factors can be 
obtained from collision rate theory and transition-state theory. Molecular 
modeling can also be used to constrain activation barriers. The remaining 
unknown kinetic parameters become adjustable parameters that are 
determined by fitting the kinetic model to experimental data. 
Experimental design is used to determine the optimal process variable 
space for kinetic parameter estimatjon. Discrimination between two 
1J mechanisms is accomplished based on how well each mechanism 
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reproduces the experimental kinetic features as process variables are 
varied. Additional experimental data may be required for discrimination 
as the optimal process variable space for mechanism discrimination may 
differ from the one used for Kinetic parameter estimation. 

(3) Kinetic parameter discrimination: Mechanisms which can no longer be 
distinguished based on experimental data (because the process variable 
space for discrimination is not accessible experimentally) can be examined 
for the reasonableness? of their values of the kinetic parameters 

(4) Sensitivity analysis! A sensitivity analysis of the overall catalytic 
performance to the|kinetic parameters of the elementary steps is carried out 
on the remaining lew reaction mechanisms which have passed steps/filters 
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(2) and (3). The most sensitive parameters are the criticaj^parameters to be 
altered 

(5) Relate critical kinetic parameter to material/properties'. This step is 
probably the most challenging as kinetiy parameters may depend on 
complex combinations of measurable ^materials properties. Postulated 
materials properties may be testedf for relevance by synthesizing and 
testing materials with varying values of that property or by calculating 
with molecular modeling theicritical kinetic parameter and the specific 
property for a series of model catalytic surfaces. The outcome of that step 
is one or several structure-property relationships that relate kinetic or 
thermodynamic parameters to a measurable materials property. 

(6) Parameter optimization: Optimum values of the kinetic parameters are 
calculated from the/microkinetic model(s) given the targeted value for the 
overall catalytic performance. These values are compared to the structure- 
property relationship to assess whether or not a materials candidate with 
the corresponding properties is plausible or not. If it is that information is 
used to select materials to be investigated in the next iteration of the TC. 

[0053] Because the reaction pathway for a caterytic process depends on the properties of the 
materials, experimental kinetic data on new materials from the TC must be re-processed 
through steps (2) through (6) to investigate possible changes in reaction pathways and to 
update the kinetic hypothesi^and the associated structure-property relationships. 
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[0054] As in microkinetic analysis, kinetic Monte-Carlo simulation can boused to simulate 
the rate of elementary surface processes. However, Monte-Cari0 simulation takes into 
account the explicit effects of surface coverage, local surface atomic composition, 
structure, and spatial arrangement. The effects of surface poisons and promoters on 
catalytic performance can also be simulated. This technique can, for example, simulate 
the essential features of the kinetics of the catalyti/hydrogenation of ethylene on Pd(100). 
Other simulations are, of course, possible. Ther method successfully describes the relative 
kinetic trends as process variables are altered and, in some cases where the metal surface 
structure is well defined, has captured some degree of quantitative accuracy as well. This 
approach is unique in relating overall catalytic performance to surface materials 
properties. While the Monte-Carlo kinetic simulation has been used to elucidate the 
effects of local structure on me catalytic process, its application to systematic virtual 
screening of catalytic surfaced is novel in the CDE. In this invention, it is used to predict 
the relative performance /activity, selectivity, life - under steady-state and/or transient 
conditions) of virtual catalytic surfaces as a function of process variables and materials 
surface properties ana to provide relative ranking of these surfaces. Surfaces that show 
promising catalytic/)roperties can then be explored further experimentally. This tool is 
aimed at guiding me scientist in the exploration and optimization of catalytic materials. 



a. 



[0055] As in the micro-kinetic analysis^the catalytic process must be broken down in a series 
of elementary steps or lumpeayelementary steps. The Monte-Carlo algorithm tracks the 

spatial and temporal change/ of all surface intermediates in order to simulate the kinetics. 
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Values for activation barriers and pre-exponential factors forthese steps are necessary 
input to the model. Pre-exponential factors can be takeri from the experimental literature 
if available or calculated from transition-state theory and collision rate theory. Activation 
barriers can be obtained from experiment; ab initio calculations or semi-empirical 
methods such as BOC theory. The latter methods, while less accurate than first-principle 
quantum chemistry calculations, are considerably faster and thus can be used for fast 
screening of catalytic surfaces and their relative ranking. The bond energies of the atomic 
species used in BOC theory are calculated from ab initio computational chemistry. 
Effects of lateral interafztions between different species as a function of distance onto the 
catalytic surface grid are separated into through-surface interactions, where the presence 
of a neighboring adsorbs species affect the local surface electronic properties, and 
through-spacer adsorbate-adsorbate interactions. The former can be quantified using BOC 
theory andfthe latter using a molecular mechanics model such as a Force-Field model for 
exampl/ 



[0056] The Monte-Carlo algorithm examines all reaction scenarios on a combination of 
different grid sites (atop, bridge, hollow) and differen.e'grid ensembles and assesses their 
probability of occurrence and rate based on site/proximity and occupancy for the reactant 
and product states of the reaction. When a scenario is possible, binding energies of 
reactants, products, and intermediates/are calculated at the appropriate surface sites. All 
possible surface scenarios are computed and stored. The overall rate of each elementary 
step is then calculated based onfthe number of possible scenarios. Each step is classified 
as being either in equilibrium or dynamic based on their rate. In the randomly generated 
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time-event, the dynamic reaction that occurs is chosen based on its relative rate as 
compared to the sum of the rates of aft dynamic processes. All equilibrium processes, 
including surface diffusion, are allowed to occur after every time. 

[0057] Validation of the simulated results against experiment for a training set of materials is 
critical because of the assumptions present in the kinetic parameter estimations. 

Machine Learning Techniques 
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[0058] Machine learning is the process of building from a set of data <yfcita structure, i.e. 
model, from which new insights, knowledge and learning can be/aerived. Two possible 
applications of machine learning in the KC of the CDE ape (1) the identification and 
ranking of critical descriptors or combination of descriptors, including chemical properties 
such as acidity, basicity, and reducibility but also^rocess variables important to the 
preparation and treatment of materials; and (2) mt selection of materials candidates or 
experimental compositional as well as prpcess parameter regions to be explored 
experimentally based on the prediction ofirthe learning model. While machine learning 
techniques have been applied earlier /or the development of heterogeneous catalytic 
materials, their success has been very limited due in part to the small amount of reliable 
data available and the choice of machine learning methods. The recent development of 
high -throughput techniques ami equipment for the synthesis and assay of catalytic 
materials is enabling the collection of much larger data sets in a much shorter time and 
thus is making the application of these data-driven techniques more important and 

[0059] There are a number of machine learning techniques available, which can be 
categorized into two classes: supervised and unsupervised learning. The selection of a 
particular technique depends on the question that the scientist is asking. Supervised 
learning allows for prediction based on the data model that is generated from the training 
set. Unsupervised learning methods generate the data model itself from the training data. 
The KC contains a toolbox of machine learning algorithms that can be selected based on 
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the inquiry of the scientist. A representative list of methods is shown in Table 2. Other 
methods could be used as well. 



Table 2. Selected Machine Learning Techniques 



♦Neural Networks 


♦Decision Trees, especially ensembles of such trees 


♦Bayesian Net 


♦Genetic Algorithms / Genetic Programming 


♦Simulated Annealing 


♦Tangled Hierarchies of Sets 


♦Recursive Partitioning 


♦Clustering 


♦Hidden Markov Models 


♦Fuzzy Methods 


♦Semantic Networks 


♦Naive Bayes 


♦Similarity Mapping 


♦Support Vector Machines 


♦Self-organizing Maps 


♦Gaussian Processes 



Knowledge Management System 

[0060] A schematic diagram of a preferred embodiment of the CDE^is shown in FIG 2. The 
CDE is based on valuable catalyst and reaction knowlea^e and an understanding that is 
used to solve catalytic problems. Easy access, retrieval, and archival of this knowledge 
improves the decision-making ability and productivity of the process. The knowledge 
system is the information retrieval system thai integrates multiple information sources and 
interfaces in a logical fashion. The ^information sources include experimental data, 
modeling, and theoretical data, frorr/literature and from proprietary corporate files. For 
example, information from kinetic modeling and molecular modeling on reaction 
elementary steps for a particular catalytic reaction is stored so that it can be re-used easily 
for a different catalytic reaction which shares one or more elementary steps with the 
previous reaction. This information may be textual, numerical, structural, or logical, and 
may reside internally o/ remotely from the database. The knowledge management system 
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preferably incudes an interactive interface for use by the researcher in each step of the 
process/ 

Description of the Catalyst Development Engine 

[0061] The role of the KC is to greatly augment the scientist's ability to unearth rapidly 
knowledge that is critical to guiding the catalyst development. These tools aim at 
maximizing learning and guide catalyst selection. Critical tasks are described further 
below. 

(1) Develop an empirical learning model 

[0062] Composition/structure, particle size, support^ftd promoters can affect the surface 
properties of catalysts. Developing trends'relating key surface properties to these 
parameters will assist the scientist in deciding what "knob" to turn in order to tune the 
catalyst properties. An analysis of existing experimental and theoretical data with its 
machine learning algorithms such as tree analysis can uncover patterns in the theoretical 
and experimental data everi in highly non-linear systems. Data for important descriptors 
are mined and correlations between performance and these descriptors are formulated. 
Descriptors include synthetic parameters as well as experimental and theoretical bulk and 
surface properties. This information can be used to construct a working hypothesis for the 
catalytic reaction, which is tested against new data and improved, at each iteration of the 
CDEoycle. 

(2) Construct a microkinetic model to screen reaction mechanisms and identify key 
kinetic paramefers 
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[0063] Microkinetic modeling will be used to narrow/fne list of plausible reaction 
mechanisms that will be investigated further by/molecular modeling. A microkinetic 
model will be built for various postulated reaction mechanisms. The performance of the 
materials will be calculated as ayfunction of temperature, reactant and product 
concentrations, for a catalyst trarning set. Mechanisms that can reproduce the general 
trends in experimental data vsnll be studied further. Sensitivity analysis will be performed 
on the rate constants o^elementary steps to identify the critical kinetic parameters and 
how they affect overall performance. These critical parameters will be related to materials 
properties using data on a series of materials in order to identify trends more easily, and 
using scientific prior knowledge and know-how. 

(3) Construct theoretical models 

[0064] Computational chemistry will be used/to provide understanding of the reaction 
mechanism and confirm/strengther^propose hypotheses (relating key catalyst properties 
with kinetic parameters of elementary steps. NDI will use ab-initio DFT periodic slab 
molecular modeling techniques/to examine the elementary surface reaction steps for the 
reaction. These methods wjlf be used to calculate binding energies, activation barriers, as 
well as rate constants fo/ elementary steps. Calculated rate constants will be inserted in 
the microkinetic modefl. Critical rate constants will be correlated with surface properties 
in order to identify key material properties. The initial calculations will be done on a 
catalyst "training^ set. The catalyst training set will be synthesized and tested to validate 
the calculationyand expose the limitations of the model. The theoretical model will be 
refined to ac^unt^oV differences between calculated and observed surface properties and 
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rate constants and to improve the model. As leaa catalysts are identified experimentally 
the modeling work will focus on these materials. 

(4) Virtual screening to gume catalyst selection 

[0065] Faster semi-empirical inethods (calibrated with ab initio calculations) are used to 
calculate kinetic parameters and insert them into a kinetic model. Catalytic performance 
will be calculated /is a function of surface composition, structure, atomic spatial 
arrangement usio/g Monte-Carlo simulation. This kinetic model can be coupled with an 
optimization Algorithm that assists in searching for the optimal catalytic surface with the 
minimum /lumber of virtual experiments. Theoretical results are validated with 
experiments. 

Example: Nitrogen Oxide Decomposition Catalyst Development 



[0066] The invention may be better understood by illustration, considering an example. In 
this example, the CDE is used to identify and optimize commercially viable NO 
decomposition catalysts for use wkn lean-burn engines. These new catalysts enable the 
rapid large-scale commerci aldose of lean-burn engines and the realization of associated 
economic and environmental benefits. Currently, catalytic converters for traditional 
combustion engines are based on three-way catalysts (TWC's), which are incompatible 
with lean-burn engines. 



[0067] The key technical barrier for a viable catalyst seems to be the higher affinity of oxygen 
for the catalytic site relative to the minority NO component - site competition. Since 
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oxygen is present in much higher concentrations, the rate of NO decomposition is too 
slow to be commercially viable. NG^ecomposition is thermodynamically favored over a 
wide temperature range of^pr^cUcal interest and the only reaction products are nitrogen 
and oxygen. The goal of the example of the present invention is aimed at identifying NO 
decomposition catalytic material that overcome the limitations of current catalysts and 
would be viable for commercial use. 

A. Review of Past Catalyst Results 
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[0068] Despite extensive research efforts over the^past 10-12 years, the currently known 
catalysts that decompose NO under the oxygen-rich conditions of the lean-burn engine 
still have severe limitations preventing their commercialization. These limitations include 



low activity, inhibition by high ox^en levels, poisoning by oxides of sulfur (SO x ), and 

;tabijitv. ' 

into perspective, the average rate of NO decomposition (for a 80% removal level) needs to 



inadequate hydrothermal stability. To put the required commercial catalyst performance 



be 1.3 x 10" 4 mol g~ l rcnn ! , calculated over the range of temperatures of operation (200- 
600°C). This rate is about approximately 50 times greater than the estimated rate on the 
most active catawst identified to date. 



[0069] Recentw studied catalysts fall into three general categories: transition metals, complex 
metal oxides, and metal ion-exchanged zeolites. The present example focuses on the 
better, but still far from adequately performing active catalytic components. Table 3 
presents an overview of some of the better performing catalysts from past work. 



Table 3. Potential NO Decomposition Catalysts Evaluated 


Catalyst 


Advantages 


Disadvantages 


Cu-ZSM-5 
(GM) 


Appears to be active for both NO 
adsorption and decomposition 


• Poor hydrothermal stability - dealumi nation 
of the zeolite 

• S0 2 poisoning 


Fe-ZSM-5 


Appears to be active for both NO 
adsorption and decomposition 


■ Poor hydrothermal stability - dealumi nation 
of the zeolite 

• S0 2 poisoning 

■ N0 2 formation 


Precious metals 


■ Good hydrothermal stability 

■ Less sensitive to S0 2 poisoning 


N0 2 formation 


Pt/alumina 


Less sensitive to S0 2 poisoning 


N0 2 formation 


Co oxide 




S0 2 poisoning 


Fluorite and 
Brownmillerite 
(Eltron Research) 


• Up to 100% NO removal 

• High concentration of O vacancies 
active for NO adsorption and 
decomposition 


Seems to rely on exhaust gas CO and HC for 
NO removal 

Work done with no 0 2 present 


erovskite 


Molecular modeling study done with 


Low NO decomposition activity 
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performance testing (Japan) 




W0 3 + Pt/alumina 
(Ford) 


• W0 3 has a positive effect on NO 
removal 

■ Test was done with S0 2 present 


• Best results were not up to desired 80% 
removal level 

■ Life information not reported 

• Reduction may significantly contribute to 
NO removal since CO and HC are present 



[0070] The problem of hydrothermal stability affects most zeolite based catalysts, especially 
Cu-ZSM-5, as set forth in M. Iwamoto and H. Hamada (Catal. Today, 10, 57 (1991)). 
The zeolite catalysts are active for NO decomposition but they are prone to dealumination 
and loss of crystallinity at high temperatures (above 973 K) in the presence of water vapor 
(typically comprising 8-10% of the exhaust). Hence, Cu-ZSM-5 is unlikely to be stable 
enough to pass the durability test of the FTP (Federal Test Procedure), as set forth in J.N. 
Armor ( Catal Today, 26, 99 (1995)). Furthermore, S0 2 poisons Cu-ZSM-5. Since the 
hydrothermal stability and poisoning problems of ZSM-5 based catalysts still exist despite 
extensive efforts to improve them, the present inventors focused initially on metal and 
metal oxide catalytic material as opposed to zeolite-based systems. However, Cu-ZSM-5 
is used as a reference catalyst in the development of our models in order to understand the 
properties of the active Cu species. 

[0071] Noble-metal catalysts are active, as evidenced by their proven performance in the 
TWCs, but they do not have sufficient/FTO decomposition activity for lean-burn engine 
application without serious modification. Supported noble-metal catalysts resist S0 2 
poisoning better than Cu-ZS^tf-5 probably because of their ability to oxidize the S0 2 to 
S0 3 . / 

B. Microkinftic Analysis 
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[0072] The observed rate of NO decomposition on Pt is proportional/o [NO]/[0 2 ]. The 
reaction is strongly inhibited by 0 2 , which competes for the adsorption sites of NO. In 
addition the presence of 0 2 in the feed leads to the oxidation of the Pt surface, and its 
reconstruction, at high temperatures. The oxidized surface is less active than Pt. The 
observed rate expression on oxidized Pt is also observed on bulk metal oxides, CuO, NiO, 
C03O4, and Fe 2 0 3 . On oxides the sites would be oxygen vacancies. 

Amirnazmi and Boudart (J. of Catalysis 39, ^83-394 (1975)) proposed the following 
mechanism: 



NO + * 
NO + NO* 
N 2 0 



NO* 

N 2 0 + O* 

N 2 + Omobiie 

o 2 + * 




(irreversible step) 

(2) 
(3) 

(4) (reversible step) 



O* + Omobile 

Assuming that step (1) is the rds, /tep (4) is quasi-equilibrated, and O* is the most abundant 
intermediate, this mechanism leads to the rate expression: 



rate = MNOMl+IQtO;,]) ~ik/K [NO]/[0 2 ] with k = k] and K = K4 



(Eq. 1) 



This mechanism invokes/the formation of a mobile oxygen atom. N 2 0 is observed as product 
of NO decomposition alt low temperatures on Pt 

Another possible medhanism for NO decomposition on metals involves the dissociation of 
adsorbed NO in adsorbed nitrogen and oxygen: 



NO + * + 



NO* 



(1) 
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NO*+* 
NO + N* 
N 2 0* 
O* + 0* 



N* + 0* 



+ N 2 0* 
-► N 2 + 0* 

0 2 + 2* 



One can derive a similar rate the following rate expression, 




PATENT 



rate = k 2 K l [NO]/(l+{K 5 [0 2 ] }°p 2 ~ k/K [N0]/[0 2 ] with k = k 2 Kj and K = K 5 (Eq. 2) 

assuming that step (2) is We rds, steps (1) and (5) are quasi-equilibrated, and O* is the most 
abundant surface specie/. There is no need to invoke mobile oxygen in this case. Based on 
this preliminary kinetic analysis, a more active catalyst would require an increase in k or a 
decrease in K. In bofh mechanisms, decreasing the heat of adsorption of 0 2 should result in a 
decrease in K. This analysis points to the oxygen binding energy as a key descriptor. 

[0073] At high temperature N 2 0 is no longer observed as a product. A simpler mechanism 
was studied in more detail with microkinetic modeling to assess whether it could explain 
experimental results collected on 5wt% Pt/alumina catalyst. The- data were collected in a 
fixed-bed reactor at atmospheric pressure, temperature between 773 and 913 K, 
concentration of NO between -200 and 1760ppm, and with 0% or 1% 0 2 in the reactor 
feed stream. Exemplary data points are illustrated in FIGs. 4a and 4b. 



The reaction mechanism includes four steps: 



NO + * 



-> NO* 



(1) Molecular adsorption of NO 
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NO*+* 
N* + N* 
nitrogen 

O* + o* 



-► N* + O* 



4 N 2 +2* 

► 

M 0 2 + 2* 



PATENT 

(2) Dissociation of adsorbed NO 

(3) Associative desorption of adsorbed 

(4) Associative desorption of adsorbed 



oxygen 
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[0074] Four parameters are required for Steps 1, 3, and 4 of the reaction scheme: (1) an 
entropy change to form the activated complex from a gas-phase species to an adsorbed 
species; (2) an enthalpy change to form the activated complex from a gas-phase species to 
an adsorbed species; (3) an entropy change to form a stable species on the surface; and (4) 
an enthalpy change to form a stable species on the surface. Parameters 1 and 2 
determine the temperature dependent rate-constant for an adsorption step and parameters 3 
and 4 determined the equilibrium constant for that given step. Step 2 requires only two 
parameters: (1) an entropy change to form the activated complex for the forward step; (2) 
an enthalpy change to form the activated complex for the forward step. The equilibrium 
constant for step 2 is determined from the gas-phase thermodynamics of the overall 
reaction to convert NO to N 2 and 0 2 [2NO □ N 2 +0 2 ]. The model utilizes 
thermodynamics properties of gas-phase reactants and products that are calculated using 
the entropies, enthalpies, and heat capacities such as those published in Yaw's (Yaws, No. 



[0075] The reaction scheme requires fourteen parameters to describe the catalytic kinetics. 
The number of parameters may be reduced further with the following assumptions. 

1. The pre-exponential factors for the adsorption steps (forward of step 1 and reverse of 
step 3 and 4) are calculated from collision rate theory assuming a sticking coefficient 
of 1.0. Accordingly, the pre-exponential factors for these steps are given by: 



35 (1999)). 



A 



1 



1.0123X10 6 
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2. The pre-exponential factor for step 2 is determined assuming no entropy change 
between the adsorbed NO species and its activated complex. 

3. Both N* and O* species are mobile species on the surface. 

4. The energetics for step 3 are obtained from ab initio molecular modeling on 
Pt(lOO): 

Activation energy for reverse step 3: E rev3 = 22 kJ/mol 
Enthalpy of reaction of step 3: AH rx 3 = 34 kJ/mol 

jfl 5. The activation barrier for forward step 1 is obtained from ab initio molecular 

n 

Rjl modeling on Pt(lOO): 

* E fOTl =21kJ/mol 

Q 

ru 
□ 

6. The activation barrier for the forward step 4 is obtained from Andersson et al. (J. 

Q 

jlj Phys. Chem. B. , Vol. 103, pl0433-10439, (1999)). 

E for3 = 21 kJ/mol 

7. Adsorbed NO* is an immobile species on the surface. 

8. The heat of 0 2 adsorption is a function of adsorbed O* coverage as proposed by 
Andersson et al: 
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Enthalpy of reaction of step 4: AH ra 4 = AH Kmcoverage (100 - ot0 o .) where 

- the heat of adsorption of 0 2 at zero coverage is obtained from molecular 
modeling on Pt(100), AH wroMvcrage = 290 kJ/mol 

- a is a scaling factor between 0 and 100 

- G 0 , is the fractional coverage of adsorbed oxygen. 

With these assumptions, four parameters are adjustable in the fitting of the experimental data. 
Their fitted values are: 

- AH = -135.0 kJ/mol 

rx, I 

- E for2 = 142.0 kJ/mol 

- a = 29.5 

- AS rx2 =-28.2J/(mol.K) 



[0076] The reaction scheme and corresponding parameters describe reasonably well the 
experimental data in the absence of oxygen in the reactor inlet (Figs. 3a and 3b). This 
model explains quite well the observed changes in the rate of NO decomposition with 
changing inlet NO concentrations. When 0 2 is added to the feed, the model over predicts 
the rate of NO decomposition (an example if which is illustrated in FIG. 4c) suggesting 
that additional data/are needed to adjust the coverage dependence of the heat of 0 2 
adsorption. 



[0077] The values of the adjustable parameters were analyzed for reasonableness. Molecular 
modeling calculations show that the heat of NO adsorption decreases with increasing 
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coverage of NO and 0 2 . The lowest value is -165 kJ/mol at combined coverage of NO 
and O equal to 0.75. However, this value does not extend to the high oxygen coverage 
predicted by the model (over 90 % coverage by adsorbed oxygen). Accordingly, the 
lower heat of adsorption calculated by the model for step 1 is reasonable. The forward 
activation energy for step 2 is higher than the 107 kJ/mol value predicted by molecular 
modeling on Pt(100), but again the latter value is valid at low coverage, and an increase in 
activation barrier is expected as coverage increases. 



(4 



[0078] A subsequent sensitivity analysis showed that the overall rate was most sensitive to the 
heat of adsorption of oxygen (step 4) followed by the activation energy for the 
dissociation of adsorbed NO (step 2). These parameters are two most important 
descriptors suggested by this analysis. More complex mechanisms need be investigated in 
a similar fashion to explain the observed formation of N 2 0 as low temperatures. 

C. Identify Key Catalyst Properties 

[0079] Based on the previous microkinetic analysis, oxygen binding energy is a key property 

affecting activity. In order to confirm this hypothesis, we attempted to correlate k/K with 

a measure of the oxygen binding energy, using published data for various catalytic 

systems. There is precedent for such an approach. The logarithm of rate constants of N 2 0 

decomposition and 0 2 exchange at a given T on many oxides correlate linearly with each 

other and with the lattice parameter of these materials for each crystal structure. Since the 

rate-determining step in the exchange reaction is the desorption of 0 2 , the rate constant for 

N 2 0 decomposition depends on the heat of adsorption of 0 2 (assuming non-activated 0 2 
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adsorption), which in turn depends on the metal-oxygen distance of the oxide. By analogy 
with these results, a correlation between k/K and the binding energy of oxygen could exist 
for NO decomposition. Using Equation. 3, values of k and K can be obtained from the plot 
of 1/r versus [0 2 ] at a given T. 



[0080] Published values of k and K for Pt/Al 2 0 3 (0.6 wt % Pt, 3.7% dispersion) and various 
oxides at 973K may be used (such as those published in Amimazmi and Boudart, J. of 
Catalysis 39, 383-394 (1975)). Values of k and K may be estimated from published 
variations of l/v t with [0 2 ] at 773K for Cu-ZSM-5 (such as Si:Al=26, Cu:Al=0.83, Cu 
content: 4.9 x 10" 4 mol g"\ as published in Y. Li and W.K. Hall J. Catal, 129, 202 
(1991)). Values of k were calculated assuming a site density of 1.19 x 10 15 cm" 2 for Pt and 
10 15 cm 2 for oxides. For the oxygen binding energy, measured by (-AH a o 2 ), one may use 

the enthalpy of formation of the oxide in the case of oxide materials (such as that 
published in CRC Hanbook of Chemistry and Physics, 74 th ed., D.R. Lide, CRC Press 
(1993)). For Pt, one may use the energy of desorption of 0 2 at high coverage because Pt 
does not form a stable oxide at those temperatures. For Cu-ZSM-5, one may used the 
enthalpy of the release of 0 2 from Cu 2 + (0 2 ") to Cu + ' , such as that calculated by the 
Density Functional Theory (DFT) quantum chemical methods (Bell et al., J. Phys. Chem. 
100, pl7582-17592 (1996)). The values of k and K are summarized in Table 4. 



l/r = (l/k[NO]) + {K/(k[NO])}[0 2 ] 



(Eq. 3) 
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Table 4. Kinetic Parameters for Various Catalysts at Two Temperatures 



T(K) 


773 


973 


Catalyst 


AH 


K 


k 


k/K 


AH 


K 


k 


k/K 




(kJ 


(atm' 1 ) 


(cm 3 


(atm cm 3 


(kJ 


(atm" 


(cm 3 


(atm cm 3 




mor 1 ) 




moP 1 s" 1 ) 


mol 1 s' 1 ) 


mol' 1 ) 


') 


mol' 1 s 1 ) 


mol' 1 s* 1 ) 


Pt/Al 2 0 3 


-201.0 


2159.4 


60200 


27.9 


-201.0 


360 


758823 


2107.8 


C03O4 


-416.2 


367.1 


414 


1.1 


-403.9 


90 


722 


8.0 


CuO 


-274.4 


496.7 


288 


0.6 


-257.5 


130 


834 


6.4 


Fe 2 0 3 










-502.7 


60 


19 


0.3 










NiO 










-401.3 


110 


30 


0.3 










Cu-ZSM-5 


-163.2 


63.5 


16396 


258.2 
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[0081] The variation of k/K with AH at 773K and 973K are shown in FIG. 5. Although the 
number of data points is limited, k/K increases as AH increases at both temperatures. 
Catalysts NiO and CuO do not fit the trend well. 



[0082] As AH increases further (or the oxygen binding energy decreases), the rate of NO 
decomposition should decrease. Since NO and 0 2 compete for adsorption on the same 
sites of these catalysts, as the affinity of these sites for oxygen decreases, their affinity for 
NO should eventually do too. Thus, an optimal value of AH likely exists and a site with 
such oxygen affinity should in principle have higher activity. Based on the data of FIG. 5, 
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materials with AH > ~-200kJ/mol should be investigated further in the next cycle of 
experiments. 

[0083] Another complementary approach for the identification of important materials 
descriptors for this reaction is the use of machine learning and pattern recognition 
techniques to unearth relationships between materials properties, preparation and process 
parameters (i.e., synthesis method, raw materials, drying temperature, etc.) and materials 
performance parameters. For example, using literature experimental values of NO 
decomposition activity and activation energy on a series of oxide materials and a list of 21 
tabulated properties of the oxides or of their respective metallic element, we applied 



independent variable with redundant information. The variables are shown in Table 5. 



|4 

Q 
Q 
01 

m 



recursive partitioning to unearth relationships between dependent and independent 



variables as well as between dependent variables. This also allowed us to identify 
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Table 5 
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Variable 


Description of Oxide MxOy 


Type j 


X1 


Atomic Number 


Discrete 


X2 


Oxide Formula 


Discrete 


X3 


Group of Element M 


Discrete 


X4 


Melting Point of Oxide (°C) 


Continuous 


X5 


Observed Activation Energy for 02 Exchange, E1 (kcal/mol) 


Continuous 


X6 


Observed Activation Energy for N20 Decomposition, ENd (kcal/mol) 


Continuous 


A / 


Hpat nf Formation of Oxirip -AH°ooo*10~ 6 (k.\/mn\ C)9) 

i leal yji rui 1 1 muui i \j\ waiuc 7 1 298 irxu/ii iui 


f^rintini \c\\ ic 

L/UI III! lUUUo 


X8 


Atomic Radius of M (Ana \ 


Continuous 


X9 


Covalent Radius of M (Ang.) | 


Continuous 


X10 


Magnetic Succeptibility (E-06 cgs) 


Continuous 


X11 


Electronegativity of Cation Mn+ 


Continuous 


X12 


First Ionization Potential of M (V) 


Continuous 


X13 


Electric Conductivity of M (E+06/ohm.cm) 


Continuous 


X14 


Thermal Conductivity of M (W/(K.cm)) 


Continuous 


X15 


Polarizability of M (E-24 cm3) 


Continuous 


X16 


Refraction Index of Oxide, Rl 


Continuous 


X17 


Polarizability of Oxide (E-24 cm3) 


Continuous 


X18 


Polarizabilityof Oxygen Anion (cm3) 


Continuous 


X19 


Ionization Enerqy of Electron in Oxide (-1/X16 4 ) 


Continuous 


X20 


Bulk Structure of Oxide 


Discrete 


Y1 


Log[Observed Rate Constant at 873K and 200 torr NO, k (cm-2.s-1)j 


Continuous 


Y2 


Observed Activation Energy for NO Decomposition, Eo (kcal/mol) 


Continuous 



[0084] The dependent variables are Yl and Y2 and the independent variables are XI through 
X21. There are no independent variables from the catalyst preparation or testing process 
in this example. The input data is shown in the Table illustrated in FIG. 6. 

Results of the recursive partitioning analysis are shown in Table 7. 
Table 7 - Average Independent Variable (Descriptor) Usage Frequency in the Construction of 
1000 Trees 
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I Y1 


Independent 


Total 


Variable 


Frequency (%) 


X5 


31.3% 


No split 


30.1% 


X20 


9.4% 


X17 


9.1% 


X3 


4.3% 


X7 


4.0% 


X8 


2.5% 


X4 


1.7% 


X15 


1.4% 


X10 


1.3% 


X1 


1.2% 


X13 


1.2% 


X16 


0.9% 


X18 


0.5% 


X19 


0.4% 


X6 


0.3% I 


X12 


0.2% 


X9 


0.1% 


X14 


0.1% 


Total 


100.0% 



Y2 


Independent 


Total 


Variable 


Frequency (%) 


X5 


61 .9% 


X6 


23.1% 


X19 


12.7% 


No Split 


0.7% 


X4 


0.5% 


X1 


0.3% 


X17 


0.3% 


X8 


0.1% 


X20 


0.1% 


X13 


0.1% 


X10 


0.1% 


X16 


0.1% 


Total 


100.0% 



[0085] In this example, the most important descriptor that describes both Yl and Y2 is clearly 
X5, the apparent activation barrier for 0 2 exchange reaction. This confirms the finding by 
Winters et al. that the fact that the ability of the material to desorb oxygen is critical. In 
the case of Yl, X5 is the only significant independent variable that describes the data; 
since the probability of no-split ranks second, all other variables are statistically 
insignificant. This is due to the small size of the dataset. In the case of Y2, X6, the 
measured apparent activation barrier for N 2 0 decomposition, and X19, the ionization 
energy of an electron in a surface oxygen ion vacancy, are significant as well. 
Decomposition of NO and decomposition of N 2 0 share common elementary steps such as 
the adsorption of NO, the desorption of 0 2 , and possibly others. Also N 2 0 is observed as 
a product of NO decomposition on Pt at low temperatures. The importance of X19 
suggests that the chemisorption of NO is connected with the transfer of an electron from 
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the material to the NO molecule. The lower X19 is, the stronger the chemisorption of NO 
and also the ability of the NO to dissociate will be. 



[0086] The fact that X7, the standard heat of formation of the oxide, which is a proxy 
parameter for the oxygen biding energy, was not significant for either Yl or Y2, prompts 
a closer examination of the data. In FIG. 7, Yl is plotted against X7. 
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[0087] The materials can be classified in two classes: In class I materials, as X7 increases, Yl 
decreases; and in class II materials, Yl seems less sensitive to X7. Class I materials 
confirms the conclusions of the microkinetic analysis and Fig. 5. In class II materials, the 
rate of NO decomposition seems to increase with increasing X7, although there is more 
scatter in the data. Class II materials are generally non-reducible oxides that have strong 
oxygen binding energy. The reaction likely operates with a different mechanism on these 
materials. This points to another region of potential interest, materials with large X7 
values, which requires further exploration. 

[0088] In this example, recursive partitioning using ensemble of trees may also be effective in 
pointing to redundant independent variables that the user may not have rationalized in the 
variable selection. For example, thermal and electrical conductivity variables (X13 and 
X14 respectively) may be highly correlated. Variables X5 and X6 may also be found to 
be closely related which correctly suggests that the decomposition of N 2 0 is also inhibited 
by adsorbed oxygen. 

D. Guide selection of materials for subsequent cycle of experiments (Virtual screening) 

45 
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[0089] In order to guide the experimental effort and minimize the number of experiments to 
be carried out, it is useful to pre-screen potential materials by calculating estimates of the 
critical kinetic and thermodynamic surface parameters (descriptors) that were identified 
previously; or more preferably by calculating estimates of the NO decomposition activity. 
This allows the scientist to look for trends in a class of virtual materials and to test the 
hypothesis put forward. Semi -empirical computational chemistry methods such as 
Extended-Huckel Theory that have been calibrated with first-principle molecular 
modeling results are useful for such parameter estimations because they are much faster 
than ab initio techniques and while they yield less accurate absolute calculation they 
generate accurate trends across materials classes. The estimated parameters can be used 
in a Monte-Carlo kinetic simulation, which takes into account kinetic elementary steps 
and the various surface sites, their concentration and their spatial arrangement. 



ry [0090] Such simulations may, for example, be carried out on Pt(100) and Rh(100) to illustrate 

Q 

jji the value and tremendous potential of this approach to catalyst discovery and development 

s 

Rl efforts. The elementary steps that would be allowed to take place on the metal surface in 

the simulation are listed in FIG. 8. The binding energies of the relevant adsorbed species, 
calculated using DFT, are shown in Table 8. Only the most probable adsorption sites are 
examined in this example. For example, nitrogen and oxygen adatoms will only occupy 
the bridge and hollow sites. Preferably, the atop sites would not be examined. The values 
for the most favorable adsorption sites are presented in italics in Table 8. The activation 
energies for the elementary steps on Pt(100) are shown in FIG. 9. For example the 
activation energy for the dissociation of NO on a bare Pt(100) surface was calculated at 
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107kJ/mol. The effects of lateral interactions (or coverage) on the energies were 
calculated using two techniques. The through-surface interactions which occur via charge 
transfer between two adsorbates through the surface. These interactions were estimated 
using Bond Order Conservation (BOC) predictions and fitting the BOC parameters to 
representative DFT results. The through space adsorbate-adsorbate interactions include 
Van der Waals forces and static electronic interactions. They were predicted using the 
Merck Molecular Force Field Model. For example, the effect of oxygen coverage on the 
oxygen binding energy is shown in FIG. 8. 
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Table 8 - Binding Energies of Adsorbed Species on Pt(100) 



Species 


Atomization 
Energy 


Atop 


Bridge 


4-hold 
Hollow 


NO 


671.5 


136.6 


214.3 


157.8 


N 






403.1 


414.4 


O 






. 394.0 


359.4 


N 2 0 


1225.3 


26.5 


11.9 




N0 2 


1073.4 


118.8 


74.4 





[0091] The simulation tracks the occurrence of elementary physicochemical steps along with 
their corresponding kinetics to follow the fate of each individual molecule that adsorbs on 
the Pt surface. The kinetics for adsorption, desorption and all surface reaction steps 
enumerated in FIG. 8 are followed as a function of time and process conditions. All sites 
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on the surface are explicitly modeled. This enables to capture the effect of the local 
composition as well as the atomic structure on the catalytic kinetics. 

[0092] The simulations start by allowing first the surface to equilibrate with respect to the gas 
phase partial pressures of the reactants. Once the initial surface structure is equilibrated 
the surface reactions are allowed to take place. Subsequently, the simulation continues for 
some time to achieve the steady state. The number of product molecules that desorb can 
be directly counted, so that an explicit measure of the turnover frequency can be 
calculated and compared to experimental values. The simulation can then be run as a 
Q virtual experiment by changing the temperature, partial pressures, surface structure, and 

cj 

01 even the metals used. Using these calculated data, the turnover frequency (TOF) for 

m 

1j steady-state formation of N 2 and N0 2 may be calculated on Pt(100) at several temperatures 

J* and [O 2 ]=60.8 torr and [NO}=0.468 torr (FIG. 11). Similar calculations may be 

Q 

j-y performed for Rh(100). The binding energies are shown in Table 9 and the TOF values 

d 

\!) are plotted in FIG. 12 for Rh(100). The Monte Carlo simulation results show that 

fU Rh(100) is inactive for NO decomposition and NO oxidation because oxygen adsorbs 

strongly to the surface and inhibits NO adsorption and reaction. 



Table 9 - Binding Energies of Adsorbed Species on Rh(100) 



Species 


Atomization 
Energy 


Atop 


Bridge 


4-hold 
Hollow 


NO 


671.5 


219.6 


263.8 


256.6 


N 




353.6 


475.4 


547.8 
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0 




376.8 


474.4 


496.7 


N 2 0 


1225.3 


53.9 






N0 2 


1073.4 


177.4 
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[0093] Catalysts of Pt and Rh supported on y-alumina may be synthesized and tested in a 
fixed-bed reactor for NO decomposition at 773 and 873K, [NO]=1.338 torr and [O 2 ]=0 
torr. The catalyst properties and corresponding TOF values for N 2 formation are shown in 
Table 10. No N 2 0 and N0 2 products were detected. The Pt catalyst was obtained from 
Alfa Aesar; the Rh catalyst was prepared by incipient wetness impregnation of a La Roche 
VGH-22 y-alumina with an aqueous solution of Rh(N0 3 ) 3 . The TOF values on Pt may be 
calculated from integrated NO conversion values assuming first order in NO, which may 
be verified experimentally. The oxygen inhibition need not be taken into account in this 
calculation. A quantitative agreement between the experimental and simulation results is 
not expected. The simulations may be carried out over model Pt(100) surfaces which 
contain distinct and active four-fold adsorption sites. The experiments may be performed 
on supported particles which likely contain a much more significant fraction of the more 
stable closed-packed surfaces. Simulation results run on the (111) surface indicate that 
this surface would be relatively inactive due to the NO dissociation barrier height which is 
greater than 200 kJ/mol. The actual measured rate over supported Pt particles is most 
likely some ensemble average from simulation results over both the Pt(100) and Pt(lll) 
surfaces. The simulations over Pt(100) and Pt(lll) surfaces therefore provide upper and 
lower bounds respectively for the actual experiments. 

[0094] Although the predicted and measured TOF values for N 2 formation do not agree, the 
Monte Carlo simulation does predict the relative decrease in activity going from Pt to Rh 
surfaces. 

Table 10 - Pt and Rh catalyst properties 
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Catalyst 


Metal Content 
(weight %) 


Metal 

Dispersion (%) 


N2 formation, TOF (s-1) 


773K 


873K 


Pt/A1203 


4.8 


29 


0.0004 


0.0015 


Rh/A1203 


0.9 


85.9 


0 


0 



[0095] These results are useful as well to the microkinetic model as it provides energetic 



values as a function of coverage. This information is very valuable to the researcher who 



can focus experiments on more promising leads as predicted by the virtual screen. Virtual 
screening saves valuable experimental resources and provides insight into the surface 
\*% chemistry. Additional surfaces of Pt and Rh and their alloys can be screened with this 

a 

O method. 

m 

fU Materials Search Strategy based on Initial Cycle ofCDE 



p [0096] The MC simulation, microkinetic analysis, and data mining results all show that both 

ru 

Q the NO activation energy and the heat of adsorption of oxygen are important in controlling 

Q the surface chemistry. To probe other materials, the researcher can use the scientific 

ru 

framework defined herein to define new material search strategies. The working 



hypothesis acts as an initial guide but will be refined as one proceeds through subsequent 



rounds of the testing and knowledge cycles. Ab initio DFT calculations, semi -empirical 



extended Huckel calculations, and literature data may be used to construct periodic trends, 



such as those illustrated in FIG. 13. FIG. 13 illustrates the challenge in finding a single 



active metal. The two parameters that appear to control the kinetics (the oxygen binding 



energy and the NO activation energy) are shown to be directly correlated. Although 



metals such as Pd and Cu bond oxygen much more weakly, they are not likely to be active 
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enough to dissociate NO. Metals such as Rh and Ru, on the other hand, can readily 
dissociate NO but they tend to readily poison. 



[0097] Based on the results from the initial run through the knowledge cycle, three new 
search strategies may be identified to circumvent the limitations found in Fig. 13. They 
include: 1) generation of active atomic structural ensembles; 2) alloy formation; 3) the 
formation of magnetic alloys. 



[0098] Initial quantum chemical results indicated that the specific atomic surface structure of 

M 

p the metal can greatly affect the nature of the transition state and hence the activation 

Q 

91 barrier for NO dissociation. FIG. 14, for example, shows the comparison of DFT and 

® 

N surface science results for the activation of NO over different ideal Pt surfaces. NO may 

hi 

be found to be inactive on the Pt(lll) surface. The surface science results would suggest 

Q 

that the more open surfaces such as Pt(210) and Pt(410) can readily activate NO. The 

Q 

U) DFT results performed in the CDE confirm that the Pt(210) and Pt(410) are more active, 

a 

Pj but the activity increase is due not only to the increase in the surface corrugation but also 

to the unique structure of the 4-fold site. In addition to the Pt(410) surface, the Pt(210) 
and Pt(100) surfaces are also fairly active. The nature of square metal atom arrangement 
seems to be unique. It provides a special site whereby NO dissociates over the 4-fold 
hollow site in the center and produces N and O products which sit at neighboring bridge 
sites. It is important to note that DFT calculations for the NO dissociation on this 4-fold 
site show that the N and O in the transition state along with the resulting N* and O* 
surface species do not share any metal atom neighbors. The transition state is such that 
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two Pt atoms stabilize N* and two other atoms stabilize O* (Fig. 15). The removal of 
metal atom sharing in both the transition states and the product states dramatically lowers 
the activation barrier. The 4-fold site is unique in that it will readily activate NO but 
should not be poisoned by the products. The previous MC simulations on Pt(100) confirm 
that effect. Indeed the calculated rate of NO decomposition on Pt(100) was much greater 
than the measured rate on supported particles. Modifiers that could stabilize this site 
would be of interest. 

[0099] Bimetallics systems is a second class of materials that could lead to increased rate of 
NO decomposition based on the CDE hypothesis. Bimetallic systems provide a mean to 

a 

flj optimize the M-O bond strength along with the activity. One method to vary the oxygen 

m 

RJ binding energy systematically is to alloy a group VIII metal with a group IB metal in 

order to form bimetallic particles. A good example is the Pd-Au system, where Au lowers 

8 

Q 

j^j the binding energy of oxygen leading to orders of magnitude increase in the rate of the K,- 

b 

y<j 0 2 reaction. Group VIII metals have sites where O2 has a higher affinity for the metal site 

Q 

ry than NO but the site is active for the catalytic decomposition of NO to N2 and 0 2 . These 

noble metals are also good oxidation catalysts and this can lead to more complicated and 
offsetting chemistry such as the oxidation of NO to NO2. Other modifiers could be used 
as well. Modifiers such as Sn and group VII elements have been used to modify the 
surface properties of group VIII metals in other reactions. There are indications in the 
literature that suggest that NO coupling can occur over PtSn alloys resulting in the 
formation of either N 2 0 (g) or N 2 (g) and O*. This provides a route in which we can alter 
the dependence of NO activation and oxygen poisoning. NO coupling provides an 
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alternative path to NO removal, which may not be as sensitive to oxygen surface coverage 
as direct NO dissociation. 



[00100] A third possible strategy is the use of magnetic alloys. The details of the 
electronic structure of the reactant and the transition state suggest that the introduction of 
a second metal with specific magnetic properties may help to lower the NO activation 
barrier. NO is unique in that it already has one electron in an antibonding NO* orbital. 
Tuning the metal surface structure so that it can inject an extra electron into the NO* 
antibonding orbital would enable a low energy path to the activation of NO and minimize 
oxygen poisoning at the same time. Indeed, the activation barrier for NO dissociation on 
Pd 3 Mn(100) alloys is lower by 60-80 kJ/mol than the value on Pd(100). Manganese 
provides the appropriate magnetic centers to offer electrons to Pd. On the other hand, Pd 
would not be very susceptible to poisoning by oxygen. 

E. Testing Cycle 

[00101] Next, a library of proposed materials is synthesized and evaluated. Scaleable, 
high-throughput experimental methods are preferable as they allow more materials to be 
synthesized and tested under a broader process condition space and a much faster rate. 
This allows for the overall CDE cycle to be shorter and for the generation of enough data 
to increase the accuracy of the data-based models. A library of materials is created by 
varying composition, raw material suppliers, synthesis methods, or synthesis process 
parameters, or others. Both the dependent variables, i.e., the activity of NO 
decomposition and product selectivity and the key descriptors (independent variables) are 
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quantified. The activity and selectivity are measured at several temperatures and several 
feed gas compositions, using steady-state or transient experiments, so that minimal kinetic 
information is available for meaningful comparison of materials. The oxygen binding 
energy is calculated from the measured energy of desorption of dioxygen by thermal 
desorption or other technique. The new data are then added to the knowledge repository 
and the CDE cycle is iterated. 



[00102] Certain portions of the invention may be performed by an automated 
processing system. Viewed externally in FIG. 16, an exemplary computer system 
designated by reference numeral 101 has a central processing unit located within a 
housing 108 and disk drives 103 and 104. Disk drives 103 and 104 are merely symbolic of 
2 a number of disk drives which might be accommodated by the computer system. 

8 Typically these would include a hard disk drive and optionally one or more floppy disk 

a 

fij drives such as 103 and/or one or more CD-ROMs, CD-Rs, CD-RWs or digital video disk 

U? (DVD) devices indicated by slot 104. The number and types of drives typically varies 

□ 

^ with different computer configurations. Disk drives 103 and 104 are in fact options, and 

they may be omitted from the computer system used in connection with the processes 
described herein. Additionally, the computer system utilized for implementing the present 
invention may be a stand-alone computer having communications capability, a computer 
connected to a network or able to communicate via a network, a handheld computing 
device, or any other form of computing device capable of carrying out equivalent 
operations. 
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[00103] 



The computer also has or is connected to or delivers signals to a display 105 



upon which graphical, video and/or alphanumeric information is displayed. The display 
may be any device capable of presenting visual images, such as a television screen, a 
computer monitor, a projection device, a handheld or other microelectronic device having 
video display capabilities, or even a device such as a headset or helmet worn by the user 
to present visual images to the user's eyes. The computer may also have or be connected 
to other means of obtaining signals to be processed. Such means of obtaining these 
signals may include any device capable of receiving images and image streams, such as 
video input and graphics cards, digital signal processing units, appropriately configured 
network connections, or any other microelectronic device having such input capabilities. 

[00104] An optional keyboard 106 and a directing device 107 such as a remote control, 
mouse, joystick, touch pad, track ball, steering wheel, remote control or any other type of 
pointing or directing device may be provided as input devices to interface with the central 
processing unit. 

[00105] FIG. 17 illustrates a block diagram of the internal hardware of the computer of 
FIG. 16. A bus 256 serves as the main information highway interconnecting the other 
components of the computer. CPU 258 is the central processing unit of the system, 
performing calculations and logic operations required to execute a program. Read only 
memory (ROM) 260 and random access memory (RAM) 262 constitute the main memory 
of the computer. 
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[00106] A disk controller 264 interfaces one or more disk drives to the system bus 256. 
These disk drives may be external or internal floppy disk drives such as 270, external or 
internal CD-ROM, CD-R, CD-RW or DVD drives such as 266, or external or internal 
hard drives 268 or other many devices. As indicated previously, these various disk drives 
and disk controllers are optional devices. 

[00107] Program instructions may be stored in the ROM 260 and/or the RAM 262. 
Optionally, program instructions may be stored on a computer readable carrier such as a 
floppy disk or a digital disk or other recording medium, a communications signal, or a 
carrier wave. 



[00108] Returning to FIG. 17, a display interface 272 permits information from the bus 
256 to be displayed on the display 248 in audio, graphic or alphanumeric format. 
Communication with external devices may optionally occur using various communication 
ports such as 274. 



[00109] In addition to the standard components of the computer, the computer also 
includes an interface 254 which allows for data input through the keyboard 250 or other 
input device and/or the directional or pointing device 252 such as a remote control, 
pointer, mouse or joystick. 



[00110] The many features and advantages of the invention are apparent from the 
detailed specification. Thus, the appended claims are intended to cover all such features 
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and advantages of the invention which fall within the true spirits and scope of the 
invention. Further, since numerous modifications and variations will readily occur to 
those skilled in the art, it is not desired to limit the invention to the exact construction and 
operation illustrated and described. Accordingly, all appropriate modifications and 
equivalents may be included within the scope of the invention. 
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